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Accurately specifying storm-time ULF wave radial diffusion in 
the radiation belts 

Stavros Dimitrakoudis, Ian R. Mann, Georgios Balasis, Constantines 
Papadimitriou,^’^ Anastasios Anastasiadis,^ loannis A. Daglis,^’^ 


Ultra-low frequency (ULE) waves can contribute to the 
transport, acceleration and loss of electrons in the radia¬ 
tion belts through inward and outward diffusion. However, 
the most appropriate parameters to use to specify the ULF 
wave diffusion rates are unknown. Empirical representations 
of diffusion coefficients often use Kp; however, specifications 
using ULE wave power offer an improved physics-based ap¬ 
proach. We use 11 years of ground-based magnetometer 
array measurements to statistically parameterise the ULE 
wave power with Kp, solar wind speed, solar wind dynamic 
pressure and Dst. We find Kp is the best single parame¬ 
ter to specify the statistical ULE wave power driving radial 
diffusion. Significantly, remarkable high energy tails exist 
in the ULF wave power distributions when expressed as a 
function of Dst. Two parameter ULF wave power specifi¬ 
cations using Dst as well as Kp provide a better statistical 
representation of storm-time radial diffusion than any single 
variable alone. 


1. Introduction 

Ultra-low frequency (ULF) waves in the magnetosphere 
have long been suggested as a likely factor affecting the ac¬ 
celeration and diffusion of electrons in the outer radiation 
belt [e.g. Fdlthammar, 1965; Schulz and Lanzerotti, 1974; 
Fei et al, 2006]. ULF wave power measurements on the 
ground or in-situ have led to various attempts to derive an 
analytic expression for their effect on the diffusion coeffi¬ 
cient of particles, as a function of some geophysical index 
or solar wind parameter [e.g. Brautigam and Albert, 2000; 
Brautigam et al., 2005; Ozeke et al., 2012, 2014]. 

The choice of which solar wind or geophysical parame¬ 
ter or parameters should be used to most accurately specify 
the ULF wave power, and thereby the radial diffusion co¬ 
efficients, is important but not yet well-understood. For 
example, the often used radial diffusion coefficient specifica¬ 
tion presented by Brautigam and Albert [2000] parameterises 
the diffusion by Kp. Similarly, according to the approach of 
Brizard and Chan [2001], the radial diffusion coefficients are 
proportional to the power spectral density of ULF waves 
(see also Fei et al. [2006]) such that parameterisations of 
ULF wave power can be used to specify the radial diffusion 
coefficients. 
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In this regard, it has been shown for example that solar 
wind speed (14™) correlates with ULF power in the Pc3 to 
5 frequency range (1.7 to 100 mHz) [e.g. Greenstadt et al, 
1979; Mathie and Mann, 2000; Pahud et al, 2009; Simms 
et al., 2010; Rae et al., 2012[, with a likely physical mech¬ 
anism being the generation of shear flow instabilities along 
the magnetopause [Cahill and Winckler, 1992; Mann et al, 
1999; Mathie and Mann, 2000]. Solar wind dynamic pres¬ 
sure changes (Pdyn) have also been found to be correlated 
with ULF wave excitation [e.g. Kepko et al., 2002; Takahashi 
and Ukhorskiy, 2007; Kessel, 2008], possibly due to the gen¬ 
eration of compressional waves that transmit their energy 
to field line resonances in the inner magnetosphere [Kivel- 
son and Southwood, 1988; Lysak and Lee, 1992; Mann et al, 
1995; Hartinger et al, 2011]. However, the nature of ULF 
wave excitation during magnetic storms is not fully under¬ 
stood. For example the morphology of the storm-time mag¬ 
netosphere, as well as the solar wind structures which drive 
storm-time ULF waves such as interplanetary coronal mass 
ejections (CMEs) or fast solar wind streams/co-rotating in¬ 
teraction regions (CIRs), may be different from the condi¬ 
tions that prevail during nominal conditions. Consequently 
the ULF wave response of the system may be different dur¬ 
ing storm and non-storm times. 

In this letter we present a study of 11 years of ground- 
based magnetometer observations of ULF waves by the IM¬ 
AGE array, binned by four parameters: Kp, Dst, solar wind 
dynamic pressure, Pdyn, and solar wind speed, Kw We 
firstly examine the distribution of the ULF wave power as 
a function of these parameters alone, to analyse which pa¬ 
rameter might be most appropriate for the statistical rep¬ 
resentation of ULF wave power. The distributions of ULF 
wave power reveal that they have a very strong high energy 
tail. We hence further analyse the ULF wave power distri¬ 
bution as a function of multiple parameters, in particular 
to examine whether the storm-time disturbance index, Dst, 
should be used in addition to other parameters, such as Kp 
or I4w which have more traditionally been used to param¬ 
eterise ULF wave power and hence radial diffusion in the 
outer radiation belt. 

In section 2 we describe our dataset and processing 
method, in section 3 we present our results, and we con¬ 
clude in section 4 with a summary and a short discussion. 

2. Data Processing 

The power spectra used in this study are derived from 
11 years of magnetic field measurements from nine stations 
in the IMAGE (International Monitor for Auroral Geomag¬ 
netic Effects) magnetometer array [Tanskanen, 2009] (1 Jan¬ 
uary 2000 to 31 December 2010): Uppsala, Nurmijarvi, 
Dombas, Oulujarvi, Rprvik, Sodankyla, Kiruna, Kevo, and 
Troms0. The particular stations analysed were chosen on 
the basis of three conditions: a) spanning a large range of 
geomagnetic latitude, b) spanning a small range of geomag¬ 
netic longitude, and c) > 95% data coverage within our 
selected time period. We analysed measurements of the D 
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(east) component of the magnetic field, sampled at 10 second 
cadence, and performed the continuous wavelet transform 
with a Morlet mother function to derive hourly power spec¬ 
tra [Balasis et ai, 2012, 2013]. The frequencies used were 
selected in the range from 0.6 to 19.85 mHz, spaced with a 
uniform linear step of 0.25 mHz (chosen for easy comparison 
with [Ozeke et al., 2012], that used fast Fourier Transform 
with a similar frequency step). These cover the Pc5 (up to 7 
mHz) and Pc4 (7 to 20 mHz) frequency ranges. To remove 
most impacts from spectral features in the Pc4-5 ULF band 
which do not correspond to waves in the equatorial plane, 
such as nightside substorm bays, we only analysed data from 
daytime hours from 0600 to 1800 MLT. We made use of the 
D-component of the ground magnetic field measurements 
because it can be mapped to the azimuthal electric field in 
space assuming an Alfvenic mode and a 90° rotation through 
the ionosphere, which is then used to calculate the electric 
term of the radial diffusion coefficient [Ozeke et al., 2012]. 
The mapping procedure we used [Ozeke et al, 2009] assumes 
a dipole geometry for the magnetic field, which is problem¬ 
atic as an approximation for the nightside of the magneto¬ 
sphere; but this is somewhat mitigated by our exclusive use 
of daytime measurements. 

The ULF wave power measurements were then binned 
with Kp, I4w, Pdyn and Dst using data from the OMNI 
database [King and Papitashvili, 2005] at hourly resolution. 
In order to robustly compare the relative efficacy of each 
binning parameter, we divided the dataset into ten equal 
deciles, ensuring that each bin for each different driving pa¬ 
rameter has the same number of data points. Table 1 shows 
the lowest and highest values of the parameters during the 
period of observation, as well as the values that denote the 
borders of their respective deciles. The total number of days 
Is 4013, since there is a five-day gap in ULF measurement 
data from early 2000, which means we have 48216 daylight 
hours of data in total, and up to 4816 hourly spectra in each 
decile. 

Since the OMNI solar wind database is not complete for 
this interval, each decile has only 4809 data points, while 
each Thyn decile has 4768 data points. The difference in bin 
sizes, with the ones for I4w being 0.14% smaller and the ones 
for Pdyn being 1% smaller than the ones for Kp and Dst, can 
be considered to be statistically negligible. The total num¬ 
ber of daylight hours for each station is further affected by 
small data gaps interspersed throughout the 11 years, with 
up to 5% fewer hours in the case of Rprvik station. Such 
gaps have no effect on the relative size of the decile bins 
for each station, and assuming the data gaps are randomly 
distributed they will only have a very small effect on the 
statistical comparisons from station to station. 


Table 1. Decile borders for each binning parameter: Kp, 
solar wind speed (I4w), solar wind dynamic pressure (Pdyn) 
and Dst. 


Decile border 

Kp 

Usw (km/s) 

Pdyn (nPa) 

Dst (riT) 

0 

0 

233 

0.03 

67 

1 

0.3 

320 

0.81 

6 

2 

0.7 

346 

1.01 

1 

3 

1 

369 

1.19 

-2 

4 

1.3 

393 

1.36 

-6 

5 

1.7 

418 

1.57 

-9 

6 

2 

448 

1.81 

-13 

7 

2.7 

485 

2.14 

-18 

8 

3 

537 

2.61 

-25 

9 

3.7 

607 

3.54 

-35 

10 

9 

1189 

79.05 

-422 


As shown below, the data demonstrate that the ULF wave 
power distributions have long high power tails, including for 
Dst. In the case of Dst we further calculate the highest ten 
percentiles, splitting the highest Dst decile further into 10 
percentile bins in order to examine the storm-time variations 
of ULF wave power. Finally, we also derive two-dimensional 
probability distributions, to both examine the utility of us¬ 
ing two parameter ULF wave specifications and to look for 
differences between ULF wave power during storm and non¬ 
storm times. 

3. Dependence of ULF wave Power on the 
Selected Parameters 

Figure 1 shows the mean PSDs and their standard errors 
in each decile for Nurmijarvi (NUR) station at L — 3.4 (ge¬ 
omagnetic longitude 102.18°), for all four parameters used 
in this study, while Figure 2 shows the same parameters for 
the Troms0 (TRO) station at L — 6.46 (geomagnetic lon¬ 
gitude 102.9°). These stations span from mid-latitude to 
close to geosynchronous orbit, with L-values which span the 
outer radiation belt. The following traits stand out (which 
also apply to the observations from the intermediate latitude 
stations and which are not shown): 

1. Binning by Kp provides the largest span in mean 
power, with smallest standard error in each decile, particu¬ 
larly in the Pc5 frequency range, and hence represents the 
best single parameter with which to specify radial diffusion. 
Binning by Usw maintains a good span and has utility for 
radial diffusion simulations as well; the deciles of Pdyn and 
Dst do not provide good discriminators of ULF wave power, 
however the ULF power in the tenth Dst decile is examined 
further below. 

2. The tenth decile of the Dst distributions is clearly sep¬ 
arated from the others, such that the mean of the 1-9 deciles 
are very similar but ULF wave power increases significantly 
once Dst < —35 nT. Pdyn shows similar behavior; indeed 
the high energy tail of the ULF wave power distribution can 
also be seen in the Kp and the I4w distributions as well. 
Perhaps more significantly, the separation of the tenth Dst 
decile from the remainder of the distribution is more pro¬ 
nounced at the lower-L station suggesting the impact of the 
storm-time dependence of ULF wave power might be rela¬ 
tively larger towards the inner edge of the outer belt. 

3. There is also evidence, especially at the high L sta¬ 
tion, for a local enhancement in the frequency spectrum of 
the ULF wave power above a simple power law, which most 
likely indicates the impact of ULF wave power accumulation 
at the local field line resonance ]cf. Rae et al., 2012]. 

This L-dependence of the ULF power is explored in more 
detail in panel (a) of Figure 3, where we have integrated the 
power spectral density profile for each decile across the Pc5 
frequency range, and derived the ratio Yi of the power car¬ 
ried by the upper decile to the total power, where i refers to 
the binning parameter. In the Pc5 frequency range, for Kp 
and Dst we see a change in Yi at around L — 4.7, and for 
Pdyn and I4w at around L = 4.3, which nevertheless have 
the same effect: a decrease of Yi with L in the outer mag¬ 
netosphere. This emphasises the strength of the high power 
tail, especially at low L. This behaviour likely reflects on 
the physical mechanisms behind the generation of the Pc4-5 
ULF waves, however, a detailed discussion of this is beyond 
the scope of the current study. Nonetheless, it shows that 
the ULF wave power distributions are dominated by a high 
power tail. 

Panel (b) of Figure 3 examines the nature of the high 
ULF wave power tail by further splitting the top decile Into 
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10 percentiles. What is incredibly clear is that not only are 
magnetic storms with large Dst magnitude associated with 
very large ULF wave power, but also that the ULF wave 
power defined on the basis of an 11 year data base remains 
well ordered by the Dst parameter even at the percent level 
in this high power tail. Overall, this suggests that studies 
of the impacts of ULF waves in the radiation belts during 
magnetic storms should additionally include the Dst index 
together with another parameter, such as Kp, when speci¬ 
fying the ULF wave diffusion rates. Such an approach may 
provide an improved ULF wave power specification for the 
purposes of modelling radial diffusion. 

Figure 4 shows occurrence probability distribution 
(panel(a)) as well as the mean ULF wave power in each bin 
as a function of pairs of deciles of two of the parameters Kp, 
Usw, Pdyn and Dst for the NUR (panel(b)) and TRO (panel 
(c)) stations. Panel (a) provides an indication of the parts 
of the distribution with reliable statistics, and in panels (b) 
and (c) the mean power is normalised to the bin with max¬ 
imum mean power. Panels (b) and (c) show that the range 
of mean ULF wave power is very large, spanning 3 orders of 
magnitude. As can be clearly seen in this Figure, Kp is an 
excellent discriminator of ULF wave power, with some evi¬ 
dence that conditions of higher I4w and Pdyn also contribute 
additionally to higher ULF wave power for given Kp. It is 
also very clear that more negative Dst also contributes to 
higher ULF power, especially at the lower L station. Recall 
that the highest value of Dst in the 10th Dst decile has the 
value of —35 nT (cf. Table 1) such that effectively all storm 
times are captured in this decile (cf. also the percent level 
discrimination of the ULF wave power as a function of Dst 
shown in panel (b) of Figure 3). Significantly, this indicates 
that the storm time magnetosphere on average has signifi¬ 
cantly elevated ULF wave power for the same activity index, 
e.g., Kp, as compared to non-storm times. 

We have also verified that the data analysis presented 
here using data from the IMAGE magnetometer array and 
the Morlet wavelet to characterise the ULF wave power re¬ 
produces the electric field diffusion coefficients reported in 
Ozeke et al. [2012] and which were derived using data largely 
from the Canadian CARISMA array [Mann et al., 2008] and 
a Fast Fourier Transform approach. The details will be re¬ 
ported in a more extensive publication elsewhere. 

4. Summary and Discussion 

We have presented a new statistical analysis of the con¬ 
nection between ULF wave activity in the magnetosphere 
and four solar wind and geomagnetic index driving parame¬ 
ters (Kp, solar wind speed, solar wind dynamic pressure and 
Dst) using ground-based magnetometer observations of ULF 
wave power spanning one solar cycle. Unlike a similar study 
undertaken by Huang et al. [2010], which evaluated the effect 
of various parameters on ULF waves in geostationary orbit, 
our focus on the D-component of multiple ground magne¬ 
tometer measurements allowed us to probe the power going 
into the electric term of radial diffusion, at L-shells covering 
the full width of the outer radiation belt. Although Huang 
et al. [2010] focused on magnetic field measurements with a 
linear binning method, which does not allow for casual com¬ 
parisons with our results even when those are for an L-shell 
corresponding to geosynchronous orbit, it is worth noting 
that they also observed a correlation of the increase of ULF 
power with increases of Kp, Dst, and Usw. Since the ULF 
wave electric field diffusion usually dominates over the mag¬ 
netic field diffusion ]e.g., Ozeke et al., 2012; Tu et al., 2012] 
then our results have significant implications for the most 
appropriate specification of the ULF wave power, and hence 
the radial diffusion rates in the outer electron radiation belt. 


We conclude that of these single parameters Kp provides the 
best capability to represent the dynamic range of ULF wave 
power, which is driven in the magnetosphere. Furthermore, 
it becomes apparent that the ULF wave power present dur¬ 
ing magnetic storms (as characterised by Dst, which mea¬ 
sures the average magnetic depression near the equator as 
a result of the ring current formed by such storms) is not 
the same as that averaged over the entire solar cycle. We 
show conclusively that the level of ULF wave power present 
in the magnetosphere is larger at storm times than at non¬ 
storm times for the same level of driving conditions, at least 
as specified by Kp (which measures the average global field 
disturbances), solar wind speed, and solar wind dynamic 
pressure drivers. As a result, we conclude that including Dst 
as an additional parameter provides an improved method for 
specifying the ULF wave power and hence the rate of ULF- 
wave radial diffusion in the radiation belts. We suggest that 
future studies should examine ULF wave power character¬ 
isations which use both Kp and Dst to specify the power 
and, as such, could be compared to studies which use the 
observed ULF wave power in individual storms to drive dif¬ 
fusion. Future studies could additionally examine the ULF 
wave power relationship to whether the magnitude of Dst is 
increasing or decreasing - thereby examining the statistical 
differences between ULF wave excitation in the main and 
recovery phases of magnetic storms. 

It has been shown in the past ]e.g. Loto’Aniu et al., 2006; 
Mann et al., 2012] that ULF wave power penetration to low 
L is connected with decreases in Dst. Here that is seen with 
greater clarity, since at low L the ULF wave power comes 
predominantly from the upper decile of Dst measurements, 
which corresponds to Dst < —35 nT, while at high L the rest 
of the deciles have a larger contribution to the overall dy¬ 
namic range of ULF wave power. Significantly, since ground 
PSD measurements of ULF waves can be directly translated 
to the electric term of the radial diffusion coefficient on the 
equatorial plane in space ]e.g. Ozeke et al., 2012], our re¬ 
sults have implications for the most appropriate parameters 
to use for characterising ULF wave driven radial diffusion. 
Accurate specification of such diffusive transport rates, es¬ 
pecially at the inner edge of the outer radiation belts, will 
likely be crucial for developing a physical understanding of 
the nature of the penetration of ULF wave power to low -L 
]cf. Loto’Aniu et al., 2006]. Such analyses may also shed 
light on the physical processes responsible for the reported 
correlations between the location of this inner boundary and 
the plasmapause ]e.g., Li et al, 2006] (itself correlated with 
Dst see, e.g., [O’Brien and Moldwin, 2003]), as well as the 
recent reports of the apparently largely impenetrable nature 
of the inner edge of the outer belt [Baker et al, 2014]. 
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Figure 1. A comparison of mean magnetic field D-component power spectral densities as measured at Nurmijarvi station 
(L=3.4), when binned by deciles, with Kp, solar wind speed, solar wind pressure, and Dst. The error bars show the 
standard error. 
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Figure 2. As in Figure 1 but for Troms0 station (L=6.46). 
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Figure 3. a) Percentage of total power in the upper decile, for the Pc5 frequency range, for all four binning parameters and 
for all stations, b) Mean magnetic field D-component power spectral densities as measured at Tromsp station (L=6.46), 
binned by Dst as in Figure 2, except that the top decile is further split into percentiles (denoted by 91 - 100). The lower 
nine deciles are plotted in grayscale (only three of them, deciles 1, 5, and 9, appear in the legend, as indicative of the 
colour progression). 
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Figure 4. a) Occurrence probability distributions for pairs of deciles of combinations of the parameters Kp, Kw, Pdyn 
and Dst. b) Mean ULF wave power for each pair of deciles, as in panel (a), for Nurmijarvi (NUR) station, c) Same as 
panel (b) but for Tromsp (TRO) station. 
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